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Abstract 

We study the effect of any uneven voltage distribution on two close cylindrical conductors with 
parallel axes that are slightly shifted in the radial and by any length in the axial direction. The 



investigation is especially motivated by certain precision measurements, such as the Satellite Test 
of the Equivalence Principle (STEP). 

> 

fN| ' By energy conservation, the force can be found as the energy gradient in the vector of the shift, 

On' 

which requires determining potential distribution and energy in the gap. The boundary value 

rn 

. problem for the potential is solved, and energy is thus found to the second order in the small 

o . 

transverse shift, and to lowest order in the gap to cylinder radius ratio. The energy consists of 
three parts: the usual capacitor part due to the uniform potential difference, the one coming from 
^ . the interaction between the voltage patches and the uniform voltage difference, and the energy of 

H : 

patch interaction, entirely independent of the uniform voltage. Patch effect forces and torques in 
the cylindrical configuration are derived and analyzed in the next two parts of this work. 
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I. INTRODUCTION 



A. Background and Motivation 

The electrostatic patch effect (PE) [1] is a nonuniform potential distribution on the surface 
of a metal. Charges in a metal move, due to a finite conductivity, until the electrostatic 
potential becomes the same at all the surface points. However, a nonuniform dipole layer 
may form on the metal surface due to impurities or microcrystal structure. If it exists, then 
the potential on the surface is no longer uniform, and the electric field is not necessarily 
perpendicular to the surface of the metal. In the (idealized) case of an isolated conductor 
the net force and torque on the body is still equal to zero, but with two metallic surfaces at 
a finite distance, the net force and/or torque on each of them does not vanish, in general. 
This effect is larger, the closer the surfaces, as first confirmed by the calculation of the patch 
effect force for two parallel conducting planes J2|. 

This calculation of the patch effect force was particularly motivated by the LISA space 
experiment to detect gravitational waves (see c.f. |3j). PE can similarly affect the accuracy 
of any other precision measurement if its set-up includes conducting surfaces in a closed 
proximity to each other. For instance, PE torques turned out one of the two major difficul- 
ties (4-6] in the analysis of data from Gravity Probe B (GP-B) Relativity Science Mission; 
the GP-B satellite was in flight in 2004-2005, to measure the relativistic drift of a gyroscope 
predicted by Einstein's general relativity (see c.f. |7j). This required theoretical calculation 
of PE torques js] for the case of two concentric spherical conductors (the lab evidence of the 
patch formation on the surfaces of the rotor and housing was found in 9|). 

The aim of this paper, consisting of three parts, is to study the force and torque due to 
voltage patches on two close cylindrical conductors. In particular, in the STEP experimental 
configuration 10h13| each test mass and its superconducting magnetic bearing form such 
a pair of cylinders. The goal of STEP is the precise (1 part in 10 18 ) measurement of the 
relative axial acceleration of a pair of coaxial test masses, so the importance of properly 
examining cylindrical patch effect is evident. 
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B. Structure of the Paper 



We determine the PE forces (and axial torque, in the final part of this paper) between 
two cylinders with parallel axes by the energy conservation argument. It implies, in general, 
that a small shift, r °, of one of the conductors relative to the other causes an electrostatic 
force between them given by (see, for instance, jhlj ]) 

where W(f°) is the electrostatic energy as a function of the shift. 

According to the formula ((!]), the force to linear order in a small shift, which is our 
ultimate goal, requires calculating W(f °) to quadratic terms in r °. Due to the specifics of 
cylindrical geometry, a more general result for an arbitrary axial shift, z°, but only small 
transverse shifts, x° = x®, y° = x% is available. So, we first find the energy in the form 

W(r ) = W (z ) + W,(z°)xl+ W^(z )xy u + O(pl), p = J(x°) 2 + (y°f . (2) 

Here and elsewhere in the paper we use the summation rule over repeated Greek indices: the 
summation always runs from 1 to 2 (meaning the directions transverse to the cylinder axes) . 
To compute the coefficients Wo(z°), W^(z°), W^i^z ), and then the force components by 
the formulas ([I]), we need the distribution of the electrostatic potential in the gap between 
the shifted cylinders to the quadratic order in the shift, 

$(f, r °) = $ (f, z°) + %(r, z°) xl + % u (r, z°) x ^ x° v + O (pg) , (3) 

since energy is a quadratic functional of the potential. 

A more rigorous view of expansions in the above formulas consists of assuming the trans- 
verse shift, po, small as compared to the nominal gap, d = b—a, between the coaxial cylinders 
(a < b are the inner and outer cylinder radii). The relevant small parameter is po/d, and 
the expansions go in dimensionless quantities x®/d, p = 1,2, so the remainder estimate 
O (po) means, in fact, O (pj/o? 3 ), etc. Also, for typical experimental conditions, such as the 
STEP configuration [ljlljjjl, the gap is much smaller than either of the radii. This justifies 
the model of infinite cylinders, and allows for a significant simplification of the results to a 
certain order in d/a. Thus two small parameters are actually involved in the problem, 

po/d < 1, d/a ~ d/b < 1 . (4) 



We work all the way to quadratic order in the first of them, and give the final answer for 
the PE forces to lowest order in the second one. However, some meaningful intermediate 
results are valid without the last or both of these assumptions. 

All said pretty much defines the structure of the paper. In the next section the boundary 
value problem (BVP) for the potential with general voltage distributions on the cylinders 
is solved, and the potential in the form ([3]) is found. Based on this, the energy represen- 
tation (jSJ) is obtained explicitly in section IIHI The details of calculations, in places rather 
complicated and cumbersome, are found in the three appendices. The derivation of PE 
forces and torques and the analysis of their properties are given in the next two parts of this 
work. The results of this part can also be used whenever the solution to the BVP for the 
Laplace equation for the domain between two cylinders is needed, such as for instance, in 
magnetostatics, thermostatics, stationary diffusion, etc. 

II. ELECTROSTATIC POTENTIAL BETWEEN TWO 
INFINITE CYLINDERS WITH PARALLEL AXES 

A. Boundary Value Problem 

We employ both Cartesian and cylindrical coordinates in two frames of the inner and outer 
cylinders as shown in figHJ In the inner, or 'primed', frame the position of a point is given 
by the vector radius r , and Cartesian coordinates {x ,y',z'} or cylindrical coordinates 
{p ', ip ', z '}. In the outer, or 'unprimed', frame the corresponding quantities are r, {x, y, z}, 
{p, if, z}. The frame origins are separated by r °, hence the primed and unprimed Cartesian 
coordinates are related by 

— * 1 — * — »0 1 1 1 / r \ 

r = r + r ; x =x + x , y =y + y , z = z + z ; (5) 

equivalenlty, in cylindrical components of the transverse shift [po is defined in (J2J)], 

x ' = x + pocosy^o, y' = y + p sin <p ; tany? = y °/x , (6) 

As alternative writing we use x = a^ , y = x 2 °, x = x\, y = x 2 , etc. 

The surfaces of the inner and outer cylinders are p' = a and p — b, respectively. They 
carry arbitrary voltage distributions, so the electrostatic potential, satisfies the Laplace 



4 



equation in the gap, 

A$ = 0, p' > a, p<b, 0<(p<2n, \z\ < oo , (7) 
and the boundary conditions of the first kind at the boundaries: 



G(<p,z), $ 



p =a 



= V- + H{ip,z); (8) 

P =b 



V~ = const (9) 



the latter are formulated using two sets of coordinates, primed and unprimed. [The potential 
distribution can also be described by using either of them, such as $(r ') = $(r*',r °), or 
$(r) = $(r, r °); the second argument emphasizes the dependence on the shift between the 
cylinders.] Boundary voltages are split in two parts: the uniform potential difference, V" 
(all the voltages in the problem are counted from the uniform voltage of the inner cylinder 
taken as zero), and the non-uniform potential distributions (patch patterns) described by 
arbitrary smooth enough functions G(<p ', z ') and H(ip, z). 

The local nature of the patch distributions is stressed by requiring 

2n oo 2tt oo 

\\G\\ 2 = J J dy'dz' \G{^',z')\ 2 < oo, \\H\\ 2 = J J d V dz\H{ V ,z)\ 2 < oo ; (10) 

— oo — oo 

these conditions are assumed valid throughout the paper. (Of course, primes can be dropped 
at the variables under the integrals, as done everywhere below). Later we will assume 
the boundary functions more smooth, expressing the additional conditions as the square 
integrability of various derivatives of G(ip , z) and H((p, z). 

The main tool in the solution of the boundary value problem (BVP) (JTj), (jSJ) is the Fourier 
transform in the axial and azimuthal variables. For any function u(ip, z) satisfying the square 
integrability condition the following representations hold: 

i oo 2tt oo 

u{ip,z) = — / dk ]T u n {k)e^ kz+ntp \ u n (k) = — / / di P dzu(i P ,z)e- i ^ + ^\ (11) 
For any two such functions u(ip, z) and v(ip, z) the useful Parceval identity holds: 

2?r 00 oo 

p p POO ^ 

{u,v) = / / d(pdzu((f,z)v*((p,z) = / dk u n{k)v* n {k) ; (12) 

J J J-oo „__„ 

-oo n ~ 00 

here and elsewhere the star denotes complex conjugation. In a particular case u = v the 
identity ([121) shows that the square of the norm, | \u\ | 2 , of a function u, defined by the formula 
ffTUj) . is equal to the square of the norm of its Fourier coefficient u n {k). 



Finally, we note that the split in the boundary conditions implies the corresponding 
representation of the potential, 



$(f) = $ u (f) + $ p (r) 



(13) 



where $" is originated by the uniform boundary voltage, and <£> p by the patch one. 



B. Solution for the Patch Potential 



According to the formulas ([7]), (JHJ) and ( TT3|) . the BVP for the patch potential is: 

A$ p = 0, p ' > a, p < b, < (f < 27T, \z\ < oo ; 



(14) 



p =a 



. , 1 roo °° / ; 

G{ V ,z) = — dk £ G n (k)e^ kz ); 

Z7T J-oo 



n=— oo 
oo 



p=6 



= — / rfA; J] tf n (fc)e 

Z7T J-oo „_ ' 



i(kz+rvp) 



(15) 
(16) 



Since the boundary functions are real, their Fourier coefficients satisfy important relations, 



G n (k) = G*_ n (~k), H n (k) = H*_ n {-k) ; 



(17) 



which are necessary and sufficient conditions for the corresponding imaginary parts to vanish. 
The dimension of functions G n (k) and H n (k) is volt ■ meter, since k is in inverse meters. 
The standard separation of variables in the Laplace equation in cylindrical coordinates 



(see c.f. [15J, Chs. 5, 6) provides the following representation of the potential: 

® P (r') = — dk \A n (k)I n (kp') + B n (k)K n (kp')} e^ kz +^ > , (18) 
2tt J-oo n= _ 00 L 

where / n (0> Kn{C) are the modified Bessel functions of the 1st and 2nd kind, respectively. 
In fact, the Macdonald's function K n (£) requires some definition for the negative values of 
its argument, which is taken here according to the parity of I n (£): the symbol K n {kp ) 
stands actually for (sign k) n K n (\k\p '). 

The unknowns A n (k), B n (k) are to be determined from the boundary conditions; the 
first of them, condition ffl5|) at the inner cylinder, can apparently be fulfilled immediately: 



A n {k)I n {ka) + B n {k)K n {ka) = G n (k), 



n = 0,±1,±2, ... . 



(19) 
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To satisfy the remaining boundary condition (fl6l) we express the ba- 
sis solutions of the Laplace equation in the primed cylindrical coordinates, 
I n (kp')e llykz +nLp \ K n (kp')e^ kz +nLp \ in terms of the set of same solutions in the unprimed 
coordinates of the outer cylinder. This re-expansion approach has been successfully used 
for solving various BVPs with the basis solutions in spherical coordinates [if]], , and 



in spherical and cylindrical coordinates 16j, 18j. The needed re-expansions of cylindrical 
solutions are readily available from 19j; the corresponding transformation formulas, along 
with all the detailed calculations, are found in Appendix A. In particular, the potential is: 

® P (r) = — dk £ \A n (k)I n (kp)+B n (k)K n (kp)]e« kz+n ^ ; (20) 

27r J -°° n=-oo L J 
00 

A n {k) = e ikz ° A m (k)I m -n(kpo)e l(m ~ n)lpo ; (21) 

m=— 00 

00 

B n (k) = e lkz ° ]T B m (k)(-ir- n I m ^n(kp )e^ m - n ^ ; (22) 

m=— 00 

n, m — 0, ±1, ±2, .... 

The boundary condition ( fl6i) at the outer cylinder can now be satisfied by requiring 

A n {k)I n {kb) + B n {k)K n {kb) = H n {k), n = 0,±l,±2,... . (23) 

Written explicitly using the expressions (l2~Pj) . (|22l) . this, along with the equation ( fl9l) . gives 
an infinite system of linear algebraic equations for the coefficients A n (k) and B n (k): 

I n {ka)A n {k) + K n {ka)B n {k) = G n {k), n = 0, ±1, ±2, . . . ; 

00 

I n {kb) Yl A m (k)I m ^ n (kp )e^ n ^ + (24) 



m=— 00 
00 



K n (A;6) £ ^(^(-^--"/^(^e^" 1 -^ = H n (k)e 



ikz° 



The solution to this system provides the patch potential for any values of the parameters 
involved. However, all we need is a perturbative 2nd-order solution in a small transverse 
shift [recall the condition (j4j) and our summation convention]: 

A n {k) = Al{k) + A»{k)xl + A^{k)xlxl + 0(p$ ; (25) 
B n {k)=Bl(k) + B»(k)xl + B^(k)xlxl + 0(pi) , n = 0,±l,±2,... . 

The system (124"]) is perfect for these perturbations: all its non-diagonal matrix elements 
are small due to the factor I m - n (kp ) ~ O (po™ - ™')- The calculations are done, and the 



expressions for all the coefficients are determined in Appendix A. They are subsequently 
simplified there to 1. o. in the second small parameter from fll]), d/a <C 1, allowing for the 
expansion of the patch potential in both primed and unprimed coordinates in the form: 



cF = $g + (xl/d) + $^ (x°Jd) (x°Jd) + O 



\iv \ jit 



(26) 



To save space, we introduce a special notation 

1 r°° 



/oo 
dk £ 
-oo 



2lT j-oo 

using it, in the inner cylinder coordinates, to 1. o. in d/a <C 1, we have: 



(27) 



G n (k) - H n (k)e- lkz ° n n (kp')e 



i (^kz 



/ i\ kz +n<p 



') = ~ t\& [G n±1 (k) - H n±1 (k)e- tkz °] n n (kp')e 



(, 1 1 

I kz +nip 



(28) 
(29) 



= --/{Vt^ (G n±2 (k) - H n±2 (k)e-^°) + 



8^/2 (G n (k) - H n (k)e 



—ikz° 



, i[kz +n V 

\l n {kp )e 



/ / 



fi,v = l,2 (30) 



[see formulas (A. 15) — (A. 17) of Appendix A]. Here 



ttn(kp') = K n (kb)I n (kp') - I n (ka)K n (kp') , 



(31) 



and terms like c ± G n ±i and c ± c ± G n ±\ should be read as the sums 



c G n ±i — c + G n+ \ + c G n -i 



c^c^ G n ± i = c + c + G n+1 + c c G n . 



(the definitions of are given in Appendix A). For \k\ — > oo, the large-argument asymp- 
totics of the Bessel functions ([3], Ch. 5) leads to 



fi„(fcp')Hr 1/2 



e -\k\(b-p ) _ e -\k\(p -a) 



(32) 



so the integrals converge and the above representations hold for a < p ' < b. 

For brevity, we omitted the factor [1 + 0(d/a)} that expresses the correction in terms 
of our small parameter d/a in all the formulas (|28p — fl3"Ul) . The estimates of the remain- 
der are only valid under some additional, as compared to f fTOj) . smoothness conditions on 
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G(ip ,z ), H((p,z), namely, the square integrability of their mixed second derivatives [see 
formulas (lAlljl . (lA12jl . Appendix A]. 

By a similar fashion, we obtain the potential in the coordinates of the outer cylinder, as 
explained at the end of Appendix A: 



a 



if n (fc) tt n (kp)e 



} 



a 



^(^) = -~ d F{[4ct (G n ± 2 (k)e ikz0 - H n±2 {k)) + 6^/2 (G n (k)e lkz ° - H n {k) 



n n {kp)e li - kz+n ^) 



p, v = 1,2 



(33) 
(34) 

(35) 



with VL n {kp) defined in (I3TI) . and the same meaning of the terms G n ±i, etc. Once again we 
dropped the factor [l+0(d/a)], which is true under the same conditions (jATTjl . (TAT2]) . Using 
the large argument asymptotics of the Bessel functions again, one shows the formulas (133]) - 
( 1351) to hold in the domain a < p < b, so that the combination of the two representations, 
in the unprimed and primed coordinates, covers the whole domain of the gap between the 
cylinders, that is, a < p , p < b. 



C. Alternative Method of Finding the Patch Potential 



The solution for the patch potential can also be obtained by the method of perturbation 



of the boundary 



20| | exploiting the first of the parameter conditions (j3J). Indeed, by the 



coordinate transformation (jSJ), using (po/p) ~ (po/ a ) "C 1, one finds 



p ' = \jp 2 + pl + 2pp Q cos(</> - if ) = p + (po/p) cos(v9 -if ) + ... . 

Hence the boundary equation p ' = a is equivalent to 

p = a- e((p) , 
e(cp) = aM %l + b^ip) xlxl + O (pj») , 



(36) 
(37) 



with the coefficients a M , b^ found explicitly. Expanding the l.h.s of the boundary condi- 
tion (fl5l) at the inner cylinder in a Taylor series according to the equality (1361) we find: 



p =a 



p=a—e 



p=a 



dp 



e + 



p=a 



1 d 2 $ p 

2~d P Y ' 



p=a 



e 2 + (pi) = G n (k) . (38) 



The potential and the boundary value function can be written in the form 



<F(f) = $g(r) + $£(r) xl + $^(f) xl x° + O (pg) 
G(<p',z') = G (ip, z) + G^, 2) + G^fa, z) x° + O (pg) 



(39) 
(40) 



Introducing the expansions (J3"7j) . (13"9"j) and (14"0"1) to the boundary condition (1551) and equating 
the coefficients at the same order on either side, we arrive at the set of the inner boundary 
conditions for the potentials of all orders: 



<T)P 
^0 



p=a 



Gn(<p, z) + 



p=a 



dp 



p, v = 1,2 ; 



p=a 



G^^.z) + 



1 d 2 $£ 

2 dp 2 



p=a 



These potentials satisfy also the following boundary conditions at the outer cylinder, as 
implied by (1TB1) and (1591: 



^0 



H(ip,z); 



P =b 



0. 



p=6 



/i, i/ = 1,2 



(41) 



and of course they are solutions to the Laplace equation in the gap a < p < b. 

Thus we have a sequence of boundary value problems that can be solved by the standard 
separation of variables in the unprimed cylindrical coordinates, obtaining the potential to the 
required order. This approach turns out much more cumbersome than that of the previous 
section. Nevertheless, it provides an important cross-check; for this reason we have carried 
it out and eventually obtained exactly the same result for the patch potential. 



D. Potential Due to Uniform Voltages and the Final Result 



It remains to determine the potential, $ u , generated by the uniform boundary voltages. 
As seen from the equations (171), (jSJ and representation ({TBI , the BVP for it is: 



A$ u = 0, p > a, p < b, < ip < 27r, \z\ < 00 



(42) 



0. 



p =a 



V 



p=b 



where V is the given constant voltage from the formula ([9]). This prob- 
lem is a two-dimensional one, since the uniform potential does not depend on 
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the axial coordinate. Its solution can be obtained by the method of sec- 
tion III Bl without even using the re-expansions of the three-dimensional cylindri- 
cal solutions. However, formally the BVP (I42p is a particular case of the BVP 
(J7|), OH]) in which the boundary data are specified through their Fourier coefficients as 



G n {k) = 0, 



H n {k) = 2nV-6{k)5 n0 



(43) 



where S(k) and 5 n o are the Dirac delta-function and the Kronecker symbol, respectively. 
With certain caution, the solution can be thus obtained from the results of section III Bl 
namely, formulas (128|) — (1311 and (13"3"j) — (1331) . The appropriate calculations, with some phys- 
ical explanations, are found in Appendix B; here we give just the resulting expressions. 
The solution of the above BVP in the gap a < p , p < b is of the form: 



= + ®l (xl/d) + $; u (xl/d) (x°Jd) + O ((po/df 
In the inner cylinder coordinates (a < p < b): 



(p /a — a/p 



p'/a 



a/p 



In the outer cylinder coordinates (a < p < b): 
r _p — a 



v 



d 



= A (c+e-*») (p/a - a/p) ; 



p — a 



(44) 

(45) 
(46) 

(47) 

(48) 
(49) 
(50) 



Note that $o(r* ) and &o(r) are identical functions of just different arguments, which, 
of course, coincide to the lowest order (p ' = p for p = 0). More importantly, the two- 
dimensional BVP fT42l) allows for a very simple exact solution in the bipolar coordinates 
{a, (3} (see e.g. [2l|, 7.3): it is a linear function of just one coordinate a, a transverse 
coordinate in the asymmetric gap. The symmetric limit (coaxial cylinders) is, however, a 
singular one for these coordinates, since the focal distance of the coordinate system tends 
to infinity. Nevertheless, it is still possible to expand the solution properly and reproduce 
exactly the above results, — another important check of our calculations. 
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III. ELECTROSTATIC ENERGY 



Denote T>l a, finite domain between the two cylinders cut at z — ±L; T)^ is then the 
whole infinite volume of the gap between the cylinders. The electrostatic energy stored in 
V L is expressed through the potential as 

W L = ^ I (V$) 2 dV = ^ / dV f(V$") 2 + 2 (V$ u • V$ p ) + (V$ P ) 2 1 , 
2 jdj, 2 Jv L L J 

according to representation (fT3"j) . The energy consists, naturally, of three terms: the one due 

to the uniform potential, the contribution coming from the interaction between the uniform 

and patch potentials, and the energy of patches. The uniform potential does not depend 

on the axial coordinate, so the uniform part of energy (and eventually the force from it, see 

section ITB"]) is proportional to the cylinder height, 2L; i.e., it is the energy per unit length 

that is finite. We thus write this part of the energy in the form: 

W U {L) = W U {L, r °) = e L f (V ± $ u ) 2 dA ; (51) 

C is the gap annulus at z = const, Vj_ is the corresponding two-dimensional gradient. 

In the two remaining contributions one can take the limit L — > oo, assuming, as usual, a 
local character of the patches: 

W int = W int^ =t [ ( V$ « . V$P ) d y . ( 52 ) 

W p = W p (f°) = ^ / (V$ p ) 2 dV . (53) 

Using the results of the previous section, we are going to calculate the quantities fl5T|) — fl53|) 
one by one in the three sections below; some details of these calculations can be found 
in Appendix C. In these calculations we systematically employ the following convenient 
formulas. 

Let T> be some domain with the boundary B, and u(r) and v (r) be some functions squarely 
integrable over T> with their derivatives. Let also u(r) satisfy the Laplace equation in T>\ 
then by integrating by parts we find 

= / vAudV= [ v^-dA - [ (Vd-Vu) dV , 
Jv jb on Jv 

where n is the outward normal to the boundary B. Thus 

/ (Vv • Vu) dV = [ v^dA, Au = in V ; (54) 
Jv Jb on 
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for v = u, in particular, 



(Vu) 2 dV = u — dA, Au = in V . 
v Jb on 



(55) 



Moreover, if both functions are harmonic in V, then we can represent the integral (154")) in 
two symmetric forms, 

/ ( Vv ■ Vu) dV = / v — dA= / u — dA, Au = Av = in £> . (56) 
Jv Jb on Jb on 

Formulas f )54|) — fl56|) hold in the space of any dimension D, in particular, for D = 2, 3. 



A. Uniform Energy 

The definition (|5T1) and formula (|55|) with w = give: 

iy u (L,f°) = e L f (V ± $ u ) 2 cM = e L I** bd^ u 
Jc Jo 

r 2n 



dp 



P =b 



e LaV I dip- 

Jo op 



[1 + 0(d/a)] 



P =b 



where we used both boundary conditions (1421) . Substituting the expansion 
potential and then carrying out the straightforward integration using expressions 
we obtain (Al is the area of the piece of cylinder with the height 2L): 

W u (L,f°) = W%{L) + W;{L) (xj/d) + W; V {L) (xl/d) (x°Jd) + O U P o/d) 



for the 

HD-dSDD, 



(57) 



W U (L) = e LaV~ / dip — ^ 
jo op 



P =b 



a / \2 
27rLe -(y-) = e 



Al(^ 



2d 



2tt <9$« 



W; U (L) = e LaV~ jj dip 



2tt /M>" 



dp 



W£(L) = eoLaV- dtp Qp 
= vrLeo-r v J = e 



(58) 



The zero order term here is, of course, the classical energy of a plane capacitor, the same as 
for the cylindrical one to l.o. in d/a; the first order term vanishes due to symmetry. 
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B. Interaction Energy 



We apply formula (|54|) to the expression (|52|) with u = $ p and v = taking into 
account both boundary conditions in the problem (|42p . which gives: 



iy mt (f u ) = e / (V$ u ■ V<F) = e aV- / cfe / d<p 



v. 



2tt ,9$P 



<9p 



[1 + 0(d/a)} 



P =b 



As in the uniform case above, we substitute here the expansion f l2"6"j) . to get: 
W mt (f°) = Wt + (x°Jd) + Wj* (arj/d) + O [(p /d) 3 



e aV~ 



oo ,2,r 

-oo jo ap 



(59) 



int 



e aV~ 



Wfi = e aV~ 



P =b 

oo 



dz / d(f—^ 

oo JO Op 



p=b 



f 2n dKv 
dz / d(p—^- 
oo jo c>p 



P =b 



The derivatives of the potential at p = 6 are found in Appendix C, formulas (IC2 
Each of them has the form of the double Fourier transform ffTTT) , 



(60) 
1. 



u(<p,z) 



I /-oo °° 

— / dk u n (k)e 



i(fcz+n</3) 



where w(<p, is any one of these derivatives. By inverting this formula and setting then 
k = and n = in the result, one obtains 

2vr oo oo 2?r 

u n (k) = — [ [ d<pdzu{<p,z)e- i{kz+nip \ f°° dz [ * d<pu{<p,z) = 2ttm (0) . 

Z7T J J J-oo JO 

-oo 

The last one is exactly the integral that stands in each of the expressions (160|) . Thus 
combining it with the mentioned formulas for the derivatives, after some not very tedious 
algebra using the relation (c~)* = c+ for the coefficients involved, we find: 



w 



int 


int 



-27T60-F- [(C o (0) - flo(0))] , 

47reo^-»[c+ (Gi(O) - #i(0)) 
a. 



(61) 



WJ* = -47reo 5 F- [c+c+ (G 2 (0) - F 2 (0))] + (G (0) - ff (0))} . 



In these calculations we have also used the property ( ITTj) as applied to the Fourier coefficient 



f n (k) = G n (k)e lkz - H n (k) ; 



(62) 



recall that 9ft(-) denotes the real part of (■). 
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C. Patch Energy 



We now apply the formula ( 155]) with u = $ p to the definition ( 153]) of the patch energy. 
Making use of the boundary conditions f[T5"]) and (j!6p . we obtain: 



W*(r°) = - / (V<F ■ V$ p ) 

2 JPoo 



dz / bdp H(p, z) — — — / dz 

-oo JO ' Op =b J-oo JO 



dz I adp G{p ,z )— — 7 
x Jo ap 



p =a 



Combining this with the expansion ( 126]) for the patch potential we arrive at the energy of 
patches written as in the formula (j2J): 



W*{f°) = Wl + W* (xj/d) + W» (xj/d) (a£/d) + O U Po /d) 



(63) 



where, to 1. o. in d/a, 



e a 



/ dz dip H(cp, z)—^- — dz dp G{p ,z 

J-oo Jo ' Op =b J-oo Jo 



dp ' 



(64) 

for £ = 0, /i, p,v, i. e., for all the functional coefficients in the expression ( 163]) . These double 
integrals can be calculated immediately employing the Parceval identity ()12p . for which one 
needs the Fourier coefficients of each of the two factors in the integrand in ( 1641 . The first 
factors H{tp, z), G(p ', z ) give no such problem, while the second ones are the derivatives of 
the potentials at the boundaries. Their expressions are found in Appendix C: formulas (102|) . 
( 1C3I) provide them at p = b, while ( 1C6I) . ( 1C7I) work for p ' = a. All these expressions have 
a structure of the double Fourier transform with the explicit Fourier coefficients. Therefore 
we can write, first for the zeroth order term in the energy expansion (163]) : 



2 



/oo °° r 
-oo „ 



/n(fc) 

d 



/n(fc) 
d 



-xkz° 



eoa 
2d 



dk £ |/„(fc)p 



(65) 



where we used the notation (162]) again, as well as the formulas (IC2I) and (1C6I) . 

Similar calculations of W£, are based on the formulas O]) . flCBl . and (IC7j) 
respectively; they exploit relations E£=-oo c~f*(k)f n ^(k) = En=-oo c~f* +1 {k)f n (k) and 
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ct resulting in 



W p 



2d 



(66) 



W p 



2d 



dk J2 [6»u/2\f n m 2 + 2M(c+4f:(k)f n+2 (k) 



Replacing now f n (k) with its expression ( 1621) in the formulas ( 1651) . ( 1661) allows us to obtain 
the final answer for the energy of the patches: 



W p 



w p 



w p 



eoa 
2d 



eoa 
' d 



eoa 
2d 



/oo 00 n 

dk £ IGnC^e*** -#„(£;) | 2 ; 
-°° n=-oo 
00 

* E S [4 ( G *ni k )e- tkZ ° - H *nW) (G n+1 (k)e^° - H n+1 (k)] 

n=—oo 

[°° dk jr, {5^/2\G n (ky kz ° -H n (k)\ 2 + 

n=-oo 

c+ct (G* n {k)e- tkz ° - K(k)) (G n+2 (ky kz ° - H n+2 (k))} } . 



(67) 



Remarkably, the expression for the zeroth order energy Wq is perfectly similar to that in 
the spherical case [see js], the last expression in the formula (3) there]. The similarity gets 
even more pronounced when rewriting the formula ( |671) for using the Parceval identity: 



W p 



eoa 
2d 



dz 



2- 



d(p [V a (ip, z) - V b ((p, z)] ; 

here V a (f, z) = G(ip, z) and Vb(<f, z) = H((p, z) are the patch voltages taken in the coordinate 
of the outer cylinder. 

The calculation of the energy is now finished; we have all its parts in exactly the form 
needed for the calculation of the electrostatic force. 
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Appendix A: Calculation of the Patch Potential 



The re-expansion formu 
coordinates are found in 



as for cylindrical solutions to the Laplace equation in shifted 



19[, 131.2.2 and 136.2.2. They read: 
I n (kp')e m * = ]T W£;po)e i(n - m)w ; 

m=— 00 

, 00 

K n (kp)e in< ? = (-l) n " m /„- m (A:po)e i(n - m) ^ K m (kp) j™" , 

m=— 00 

where po and <po are the polar components of the shift as defined in (j2J) and @. Therefore 

00 , 
X; [A n (k)I n (kp') + B n (k)K n (kp')]e m * = 

n=— 00 
00 00 

X I n (kp)e in *> Im-n(kp )e^ m - n ^°A m (k) + 

n=~oo ■m=—oo 

00 00 

X K n {kp)e in * X (-l) m -"/ m _ n (A;po)e i(m - n)vo S w (A ; ), 

n=— 00 m=— 00 

where we have changed the order of summation and then switched the indeces m and n. 
When substituted in the representation ffl8l) of the patch potential $ p in primed coordinates 
taking into account z = z + z°, this gives exactly the representation ( 1201) — (1221) in the 
unprimed coordinates of the outer cylinder. 

We are now to construct the 2nd order perturbative solution (1251) . 

*X + O0*), (Al) 



-A n (ky 






+ 


-A^k)- 






.B n {k)_ 




-B° n (k)_ 




-B%(k)_ 




-B$?(k). 



to the infinite system of linear algebraic equations fl2M . As mentioned in section III Bl 
its matrix coefficients contain the factor I m -n(kpo) ~ O(p m ~"'), so only the terms 
m = n, n ± 1, n ± 2 should be taken into account, to the 2nd order. Using thus 
the formulas (all the needed results on Bessel functions can be found in 15J, Ch. 5) 

I (kp ) = 1 + (kp ) 2 /4 + O ((A;p ) 4 ) , i±i(*Po) = (*Po)/2 + O ((kp f 

I± 2 (kp ) = (kp ) 2 /8 + O((kp )*) , 

as well as the convenient representations (cf =0.5, 62 = 0.5i) 



„,o. 



(po/2)e 



(p /2)V^ = <c+^^, 



we find: 



Jo(fcpo) = 1 + -rS^xl xl + O ((kpof) ; 



(A2) 



(A3) 
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h(k Po )e^ = kc+xl + O ((k Po f) , U(kp )e-^ = fcc^xj + O ((kp ) 3 ) , c~ = (c+)' 



h(kpo)e 



2itpo 



2 C ^ 



xJx°+0((fc Po ) 4 ) , I- 2 (^o)e- 



-2l(/30 



Here we introduced the notation c~ for a uniform writing. Substituting expressions 
and ( 1A3|) in the system (|2"%|) allows us to obtain the following sequence of equations for the 
unknown coefficients in expansions flAll) : 

I n (ka)Ai(k) +K n (ka)Bi(k) = G*{k) , 
I n (kb)Ai(k) +K n (kb)Bi(k) = Hi(k), 
where n — 0, ±1, ±2, . . ., and the index £ assumes three values £ = 0, /J, /zz/, with /x, i/ = 1, 2. 
The right hand sides here are: 



G° n (k) = G n (k), 



H°{k) = H n {k) e 



-ikz u 



-k c 



± 



<%{k) = GT = ; 

I„(*6) i^iC*) B0 ±1 (fc) 



(A4) 



2 



= - k [i n (kb)c±A: ±1 (k)-K n (kb)c±B: ±1 (k) 

5, 



Ukb) [c^ ct A» n±l (k) + ^ A» n (k) +K n (kb) Ic* ^ B» ±1 (fc) + 



and terms like c ± A n ±i and c ± c ± Ai±i should be read as 

t^AiiiW = c + A n+1 (A;) + cTA n _ 1 (A;), t^c* A»±i(A0 = c + c + A n+ i(fc) + c _ c _ A n ^(k) 

Solving the above linear systems leads to the following set of answers (£ = p, pv): 

G n (k)K n (kb) - H n (k)K n (ka)e^ kz " 



A° n (k) = 
B° n (k) 



Ai(k) 



Ht(k)K n (kc 



D n (k) 



D n (k) 

H n {k)I n {ka)e^ kz ° - G n (k)I n (kb) 
D n (k) 
B ^ = Ht(k)I n (ka) 



D n (k) 



(AS) 
(A6) 



where 



D n (k) = K n (kb) I n (ka) - K n (ka) I n (kb) . (A7) 

Expressions (1A5I) — HA7[) and (1 A 1 1) provide the solution for the patch potential in the form 
([18]) to the 2nd order in po/d and any value of d/a. We now simplify all the expressions to 
1. o. in this small second parameter. To do this, we employ the Taylor expansions: 

\2" 



I n (kb) = I n (ka) + I n (ka) kd + O (kd) 
K n (kb) = K n (ka) + K{ka) kd + O \(kd) 



a + d, d — > ; 



(A8) 
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where the primes denote the derivatives with respect to the whole argument. With the help 
of these expansions and the known formula for the Wronskian of the Bessel functions, 



w(i n (z),K n (t)) = -i/z, 

we simplify the denominator ( 1A7|) to lowest order in d/a as: 



(A9) 



D n {k) = kd I n (ka)K n (ka) — K n (ka)I n (ka) 



d 

+ - 

a 



d 
a 



l + o(-j +0 ([n(ka)~ 



(A10) 



The second term in the remainder could be dropped if k and n were bounded. But in 
our formulas both run from minus to plus infinity, so we need to keep this remainder to 
make the asymptotic expressions uniform for the whole range of k and n. However, when 
substituted in the sums over n and integrals over k, as in formula ffl8|) . these terms give rise 
to contributions 0((d/a) 2 ) provided that the proper sums/integrals converge, 



/OO °° POO 00 
dk k 2 n 2 \G n (k)\ 2 < 00, / dk k 2 n 2 \H n {k)\ 2 < 
-00 J— 00 



00 . 



By the Parceval identity (fl2l) . the latter conditions are equivalent to 

2n 00 2ir 00 

J J d^p'dz' \d 2 G/dif'dz')\ 2 < 00, J J difdz\d 2 H/difdz)\ 2 < 00 



(All) 



(A12) 



-00 



-00 



Moreover, in the expressions of A° n1 A%, A£ u we use the expansion 

K n (ka) = K n (kb) - K' n {kd) kd + O [(kd) ' 

while for B® we employ the first of the formulas flA8j) . This is done, in fact, to provide the 
proper convergence of integrals over k and series in n at infinity, see section III BL 
For the zero order coefficients, we thus have: 

2N 



A° n (k) = -- { K n (kb) G n (k) - H n (k)e 



ikz° 



(t\ +0 \[n{ka)- a 



B° n (k) = ~ I I n (ka) [G n (k) - H n (k)e-* kz °] + O (~) + O ( U(ka)~ 



(A13) 



From now on, all our asymptotic formulas have the same two-term remainder estimates, 
which we simply drop. 
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By the formulas (IA13I) and definition (1A4[) . we obtain thus 

H n(k) = -A \G n ±i(k) - H n±1 (k)e- lkz ° 



d 



and the expressions (1A6I) become 



A^k) = -^K n (kb) cj [G n±1 (k) - H n±1 (k)e- lkz ° 
K(k) = [G„±i(A;) - H n±1 (k)e- lkz " 



(A14) 



By the same token, formula (1A4I) simplifies to 



~G n ± 2 (k)-H n±2 (k)e- lkz ° 
and the coefficients (1A6j) are written as 



+ 



G n {k)-H n {k)e 



kz ' 



Kik) 



a 
"d 



, ,K n {kb) 



B^{k) 



^\ln(ka) 



A 4 V 



- H n±2 (k)e 



G n±2 {k) - H n±2 (k)e 



-i kz° 



-ikz° 



6 



Y 



G n {k) - H n {k)e 



ikz° 



(A15) 



6 



ft V 

T 



G n (k) - H n (k)e 



-ikz° 



Expressions f)A13[) — (IA15[) . when substituted in the representation (fT8"j) . provide exactly the 
formulas ( 128]) - (1301 for the potential in the primed coordinates. 

To obtain the potential in the unprimed coordinates we need, as seen from the formula 
( I20|) . to compute the coefficients A n , B n to the 2nd order in the transverse shift, in a complete 
similarity with the formula (I A 1 [) : 



A n {k) 
B n (k) 



A° n (k) 
B° n (k) 



+ 



K{k) 

BH(k) 



K v {k) 
B£>(k) 



4 4 + o (pI) 



(A16) 



To get them, we simply use the definitions (l2Tj) and ( 122|) of the coefficients A n , B n , and 
introduce there the expansions (1A1[) with the known coefficients (1A13|) — ( IA15p . as well as 
the expansions (1A3|) . This is a tedious and cumbersome, but straightforward calculation, 
which ends with the representations (j33|) - (l35l) . 
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Appendix B: Calculation of the Potential due to 
the Uniform Boundary Voltages 



The uniform potential is obtained substituting the Fourier coefficients of the boundary 
functions (JI5j) . 

G n (k) = , H n (k) = 2nV-6(k)6 n0 , 

in the formulas found for the patch potential, (I28l) -( l30l) and (|33l) -( l35l) . The zeroth order 
(coaxial configuration) potential in the primed coordinates then is [recall the definitions ( 13T|) 
and (T27j)]: 

*o(0 = \K (kb)I (kp') - I (ka)K (kp') . 

Since Io(£) — >■ 1, ^o(0 ~ l n (2/0 when £ — )■ 0, we furthermore obtain 



$«(f ') = § y- ln(p =V- [ln(p '/a) + 0(d/a) 
a a L 



(Bl) 



The logarithm above is a remnant of the exact solution for the coaxial case, 
[V~ \n(p' /a)/ ln(b/a)], with the logarithm in the denominator taken to lowest order in 
d/a « 1 (recall that d = b — a). Moreover, the upper log may be simplified in a simi- 
lar way: 



ln{p' /a) = In 1 



p — a 



p — a 



[1 + 0{d/a)} 



a l a 

the potential in a narrow capacitor does not feel the curvature of the electrodes to 1. o., 
therefore it is a linear function of the transverse coordinate. Introducing the last expression 
to the formula fIBip we finally obtain: 



') = V- [1 + 0(d/a)} . 



(B2) 



By the same token, the linear and quadratic potentials are obtained from expressions 
29D and (El as: 



K^')--l V ~ l ^M^ ltP ' [ K Ti{kb)lTi{kp')-I Tl {ka)K Tl {kp') 



and 



(ff 



-V- lim c^ 2 ^ 



K T2 (kb)I T2 (kp ) - I T2 (ka)K T2 (kp )J + 
6^/2 [K (kb)I (kp') - I (ka)K (kp')}} . 
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Since I_ n (0 = ^n(£)j ^-n(0 = these may be written as (3?(-) is the real part of (•)): 

(B3) 



, 7" lim ^ 23fJ(c+e-^ ) Kx(kb)h(kp ) - h(ka)KAkp' 



^(f') = ^-Hmj23ff( C Jc+e-^ 



K 2 (kb)I 2 (kp ) - I 2 (ka)K 2 (kp' 



+ 



<W2 [i^Ko^P ') - I (ka)K (kp ')] } . (B4) 
The near zero asymptotics of Bessel functions allows one to transform (IB3I) . flB4j) to: 



*»(T'') = --V-{»( C +e-^ )[p76-o/p' 



~V- {sR(c+e-^') [p'/o-tt/p'l +0(d/a) 



(B5) 



(p'/b) a -(a/p')* 



2d 1 M 



(p7 a ) 2 - ( a /p') 2 



5 



+ 5^ In (p 

p' — a 



/IV 



+ O (d/a) 



(B6) 



where we used the same transformation of a logarithm as when deriving formula ( 1B2j) from 
expression ( 1B1I) . 

Formulas (IB2|) . (IB5|) . and (IB6P are the final expressions for the uniform potential in the 
inner cylinder coordinates as given in formulas (|45p - (|47j) . The potential in the unprimed 
coordinates is now gotten in exactly the same way without any new difficulties: one just 
starts, respectively, from the expressions (l33]) - (l35l) of the patch potential in the unprimed 
coordinates, and follows the steps described above; in this way, the final results (l48l)-(!50l) 
are found. 



Appendix C: Calculation of the Energy 



Here we provide the intermediate results needed for the energy computation. The first of 
them is the derivatives d^/dp, d^/dp, d<& v /dp at the outer cylinder boundary p = b re- 
quired to calculate W tnt by the formulas (|60|) . First of all we differentiate in p the expressions 
fl33|) -fl35l). i.e., the function Q n (kp) in them. The result is: 



^h^Q = ^ [ Kn (kb)I n (k£) - I n (ka)K n (k£)] = k 



AT n (fc6)/ n '(fcO-/n(fca)<(fcO , (CI) 
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where the prime indicates the derivative with respect to the whole argument, as usual. At 
the outer boundary, exploiting the relation 

I n (ka) = I n (kb) - kdl' n (ka) + O ((kd) 2 ) , 

and the formula (1A9j) for the Wronskian, we can write the derivative (ICip 



dn n (kO 



-W(I n (k£), K n {ki))\K=k b + 0[(kdf 



\ a , 



Combining this with the representations (1551) — (|35p . we obtain: 



dp 



P =b 



--T{fn{k) e« kz+n ^}; 



d<$> p , 



dp 



P =b 



dp 



P =b 



- d F{[4 C »fn±2{k) + (W 2 ) f n (k)} e« kz+ ^} . 



(C2) 
(C3) 



Here J 7 is the integration — summation operator defined in ( 1271) . and, for brevity, 

f n (k) = G n (k)e lkz ° - H n (k) . 



(C4) 



The correction in the formulas (1C2|) . (IC3I) is given by [1 + 0(d/a)} under the conditions 

dk \k\ 3 \G n (k)\ < oo, dk \k\ 3 \H n (k)\ < oo . (C5) 

-oo "=-oo _ QO n=-oo 

The derivatives <9$g/<9p', <9$£/<9p', dp ' at the inner cylinder p = a are also 

needed, by the expressions ( )64|) . They are derived from the formulas (128]) — ( !30l) exactly as 
above, the only significant variation is in transforming the formula (Id p . Using the second 
of expansions (1A8|) and formula 0A9p . we find: 

i r /d^ 2 

-W(I n (k£),K n (kZ))\kt=ka + 0[(kd) 2 ? ' >/( 



dn n (kO 



dz 



k 



+ 



3„2 



k A a 



With this small change, the desired formulas are found as: 



dp ' 



= -J 7 

p =a 

,9$? (f' 

d 



f n (k) e »(*(* '-«")+«*") 



dp ' 



J" 



<9p ' 



-7 



p =a 

ctc±f n±2 (k) + (5,j2)f n (k) 



z f n±1 (k) e <fc(*'-*°)+^') 



(C6) 
(C7) 



p =a 



with the same notations and same order of the remainder as in the formulas (jC2p . ( 1C3 
assuming again condition (1C5h . 
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FIG. 1: Geometry of the problem and coordinate systems 
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